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POSSIBLE STATISTICS OF TWO COUPLED RANDOM FIELDS: 
APPLICATION TO PASSIVE SCALAR 


B. DUBRULLE* AND GUOWEI HE+ 

Abstract. We use the relativity postulate of scale invariance to derive the similarity transformations 
between two coupled scale-invariant random fields at different scales. We find the equations leading to 
the scaling exponents. This formulation is applied to the case of passive scalars advected i) by a random 
Gaussian velocity field; and ii) by a turbulent velocity field. In the Gaussian case, we show that the passive 
scalar increments follow a log-Levy distribution generalizing Kraichnan’s solution and, in an appropriate 
limit, a log-normal distribution. In the turbulent case, we show that when the velocity increments follow 
a log-Poisson statistics, the passive scalar increments follow a statistics close to log-Poisson. This result 
explains the experimental observations of Ruiz et a 1. about the temperature increments. 

Key words, scale invariance, scale transformation, scaling, turbulent scalar field 

Subject classification. Fluid Mechanics 

1. Introduction. Scale invariance [1] usually refers to physical fields keeping the same properties at 
different scales. Therefore, scale invariance is often associated with the existence of power scaling laws. In 
this case, the central characteristic quantities of the system are the scaling exponents. When a random field 
is statistically scale invariant, each of its structure functions follows a power law, and there is a whole family 
of exponents associated with its statistics. 

To understand the underlying properties associated with the scaling exponents one can follow two paths: 
first, try to measure or compute the scaling exponents with as high a precision as possible. This is the path 
usually followed in critical phenomena, where scaling exponents are few, and appear to be universal, due to 
some underlying conformal symmetry [2]. In turbulence, this path seems more difficult to follow, because the 
number of scaling exponents to measure or compute, in order to know the statistics of the field, is infinite. 
Therefore, parallel to the measurement or actual computation path, there seem to be the necessity of a more 
formal path, devoted to understand link between the various scaling exponents. That way, the measure 
or computation of only a few scaling exponents would hopefully enable the unraveling of the statistical 
nature of the system. A natural way to establish link between scaling exponents is to use the symmetry 
responsible for the existence of the scaling exponents, namely the scale invariance. This was the basis of 
conformal theory [2] or multifractal theory [3]. More recently, Dubrulle et a 1. [1, 4, 5, 6] found that new 
constraints on scaling exponents stemming from scale invariance could be found via an analogy between 
scale invariance and relativistic dynamics inspired by a theory of Nottale [7] and of Pocheau [8]. This can 
be used to classify the possible forms of scaling exponents according to the fixed points of the similarity 
transformation. However, these results were restricted to one random scale invariant field (ID case). In 
many physical situations, the dynamics involves at least two coupled random fields, such as transports of 
passive scalars (coupling of velocity and scalar), magnetohydrodynamics (coupling of velocity and magnetic 
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field) and turbulence (coupling of transverse and longitudinal velocity increments). The goal of the present 
paper is therefore to extend some of these results to the case of two coupled scale invariant random fields 
(2D case) and use the 2D results to explore scaling of passive scalars. 

The transport of passive scalars, for example the temperature 0 , is governed by the diffusion equation: 

(1.1) d t 0 + (u • Y)$ = kA9 , 

where the solenoidal velocity field u is either taken to be random or governed by the Navier-Stokes equation, 
k is molecular diffusivity. We are interested in scaling of passive scalar structure functions for given scaling 
of velocity structure functions: 

S e n (£) = ((8(x + £)- e(x)) n ) oc 

(1.2) S"(*) = <(u(s + £)- u(x)) n ) oc f (n \ 

where ((n) and £(n) are the scaling exponents of velocity structure functions and scalar structure functions 
respectively. We will use 2D similarity transformations to derive these possible scaling exponents. 

This paper is organized as follows. After useful preliminaries regarding scale invariance (Section 2), the 
2D similarity transformation of scale invariant fields is presented in Section 3. The results are applied to the 
case of the transport of a passive scalar by a Gaussian velocity field in section 4 and by a turbulent velocity 
field in section 5. We conclude the paper in section 6. 

2. Preliminaries on scale invariance. Scale invariance is a commonly used concept but sometimes, 
rigorous definitions are lacking. In this section, we give a brief review on the rigorous definition of scale 
dilation (or contraction), global and local scale invariance. Then, we discuss the existence of similarity 
transformation between measurements of physical fields at different scales due to scale invariance. 

Let (tp((x),ipt(x)) be a two-dimensional positive homogeneous random field [9] and dependent on scale 
l. We define a family of scale dilation Sh,h' (A) via: 

(2.1) Sh,h'( A) : l — > f \ h (p(,\ h ipe- 

A two dimensional homogeneous random field is said to be globally scale invariant or, simply, scale invariant, 
if (y>f(x), Vv(x)) and (A) (^ (x), ipe(x)) have the same statistical properties. 

In log- variables, T = ln£, X = In(^), Y = ln(^), the globally scale invariance amounts to a transla- 
tional invariance. This invariance is connected with the possibility to multiply any characteristic quantity 
by an arbitrary constant, i.e. with the possibility to perform arbitrary changes of units in the system [1], It 
is a global scale invariance that leads to the existence of power laws of the moments of random homogeneous 
fields: 

(2.2) (v>r)«^ (B) , (</’D oc £ c(n) . 

Scale invariance can be also considered as local [8]. Local scale invariance is connected with the invariance 
of the system with respect to resolution changes. Indeed, consider a discrete slicing of the scale space l < L 
in the form: 


(2.3) 


= L = , %, 

<P% = A Vo, i>i = A Vo, 
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where, , (or equiv alently N) and A characterize the resolution of the observations, tpo and tf’o 316 the units of 
the observations. By local scale invariance, we require that the physical properties of the system be invariant 
under changes of resolution, i.e. under changes of N. In fact, changes of N can be achieved by dilation of 
,: which transform * into ia. It is also a transformation which transforms £i/£j into (£i/£j) a , and 

< j>i/<l>j into (<f>i/<j>j) a - Therefore: 

A local scale dilation can be seen as a transformation for arbitrary a: 

(2-4) £/£ 0 — t (l/£o) a -,(p/(po,4’/4’o {vlVo) a , (4’/4’o) a , 

A two dimensional homogeneous random field is said to be locally scale invariant if the physical fields 
are statistically invariant under local scale dilation. 


Local scale invariance can be also seen as a change of rationalized unit system because it amounts to 
multiply to, f>o and tpo by £o{£l£o) l ~ a , t Po( t p/ t Po) 1 ~ a and V’o(V’/V’o) 1_a - Therefore, global scale invariance 
implies local scale invariance. 

When the field (y>f(x),r/v(x)) is globally scale invariant, all these changes of units can be characterized 
by only one number: the exponent of their power law. This number labels the system of units and varies with 
it. It is not an absolute quantity. All these remarks have been so far considered as rather obvious because 
they do not lead to any specially useful information when applied to deterministic fields. For random fields, 
however, < (p a > and < <j> a > can be very different, and the application of our remarks can lead to very 
interesting consequences about the exponents. In fact, they are not necessary additive anymore and can 
obey more complicated composition law, which can be determined only from the requirement of local scale 
invariance, as was shown in [6]. We now define the formalism which can be used to derive these composition 
laws in the case of two coupled fields. This is done by a. natural generalization of [1, 6] which we refer to for 
more details. 


3. Formalism and derivation of the composition laws. 

3.1. Formalism. As noted previously, scale invariance amounts to a translational invariance when 
logarithm of scale and fields are considered. We also would like to introduce a field variable which is 
deterministic but keeps track of the possible different realizations of the process. A good way to reconcile 
these two requirements is to introduce the log- variables [6]: 


(3.1) 


ln(l/4) 

ln(FQ 


logic (^o), 


( ln(<p ( /0 o ) (4>(/<f>o) n ) 
o) n ) 

( in (i’f/ii’o) (j'di'oY ) 
((4’e/4’o) p ) 


where K is the basis of the logarithm, £ 0 is an arbitrary length unit, </> 0 , tpo two arbitrary field units (pos- 
sibly random and/or scale dependent [1]), and n and p are real numbers. We note that ln(Af), p and n 
labels the system of subunits used to measure the length and fields. Then, globally scale invariance implies 
homogeneity in these log-coordinates. Local scale invariance implies the existence of a group transformation 
connecting different measures of the field at different scales. Formally, this means that we can find a simi- 
larity transformation between (log K (£/£o)i A'„, Y p ) and (\og K , (£/ £q), X n i ,Y p i). Since £ is deterministic, it is 
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equivalent to change K or l to change T — ln/ t -(£), so we can as well fix K e.g. K = e and let £ vary. The 
transformation is now between (T = \n(£/£ 0 ),X n ,Y p ) and (T 1 = In (£' /£' 0 ),X n i ,Y p i). We call it the similarity 
transformation. 

The derivatives U(n) = dX n /dT,V(p) = dY p /dT are analogous to the relative velocities. In our formal- 
ism, they are just the local multifractal exponents [6]. 

3.2. Two-dimensional similarity transformation. The possible similarity transformation can be 
derived provided the system is both globally and locally scale invariant. Global scale invariance implies 
that the transformation are linear [1], Since nothing forbids a coupling between fields and scales, their most 
general representation is a 3 x 3 matrix coupling X, Y and T. By local scale invariance, the set of such 
3x3 matrices must form a (at least semi) group, depending on two parameters U and V [1]. Without loss 
of generality, they may be written as: 


(3.2) 


A B , U 
D E , V , 
GU HV , 


where A , B , D, E, G, H , U, V and , are some unkno wn functions of U and V to be determined. According 
to the transformation (3.2), a field characterized by A' = Y = 0 for any T in one reference field R = (n, p) 
is transformed into X' = UT', Y' = VT' in another reference field R' = (n 1 ,p'). This shows that (U,V) are 
characteristic numbers labeling the change from (n,p) to (n',p'). It is then natural to adopt U and V as the 
group parameters. 

To completely determine A,...,, as a function of (U,V), an additional assumption is needed. From 
physical considerations, we choose to postulate that the group of similarity transformations is commutative. 
This means that it is equivalent to perform a. transformation at scale £ then £' or £' then t. This property is 
generic in the group of scale dilation S\(h , ft'), at any given h and ft', which motivates our choice. Moreover, 
it is a natural extension of the commutativity property in ID [1], Note that in such case, the commutativity 
is guaranteed by the ID Lie group structure and need not be postulated. Note also that in special relativity, 
the third additional postulate is isotropy. This leads to a representation in term of rotations, which is NOT 
commutative in D > 1. There is therefore no hope to find the Lorentz group as a special case of our similarity 
transformations, contrary to the ID case. 


With this additional postulate, the form of the matrix (3.2) can be completely specified. The proof is 
given in Appendix A. It is: 

A = (l-aU-jV), B = -(jU + dV), 

(3.3) D = —c(jU + dV ), E = (l-cdU- bV), 

G == -ft, H = -k, 

where , is a function of U and V obeying the transformation rule: 

(3.4) ,( U", V") = ,( U',V'),( U,V){ 1 - hUU' - kVV'), 

with U", U',...V following (3.7). We were not able to find the exact expression of , in the general case, 
but in the simple cases we study in the next two sections (where ft = ft = 0), , = 1 is alw ays a solution. In 
general, we postulate that , is a quadratic form in U and V, but the explicit value of the coefficient is not 
needed in the derivation of the composition law. 
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The coefficients a, b , c, d, h , fc and j are parameters depending on the physical model. They satisfy 


(3.5) h = ck, — k + ad + jb = j 2 + d 2 c 

They are real if the scale symmetry is continuous [1]. 

Composing two transformations from reference fields R = (n,p) to R' = ( n',p ') and from R' = ( n',p ') 
to R" = ( n",p "), we obtain: 

(3.6) (Ur\R" ,Vr\R") = (Ur\ri , Vr\ R! ) ® (Uri\r», Vr'\R") 
where <g> is the composition law: 

U + U' - aUU' - dVV - j(U'V + UV') 

1 - hUU' - kVV' 

V + V' - cjUU' - bVV - cd(U'V + UV') 

1 - hUU' - kVV' 

Like in ID, the composition law can be characterized by its fixed points, which depend on the parameters 
a,b,c,d,h,k,j. It can also be used to determine the scaling exponents, by an appropriate choice of R, R 1 
and R" . In both cases, the general computation is rather cumbersome and complicated by the large number 
of degenerate situations (collapsing or vanishing fixed points). We treat in the following sections two special 
cases where the computations are tractable. We emphasize that the case c = 0 implies no influence of A' 
over Y and is relevant to the case of the passive scalar dynamics when X is taken as the passive scalar 
log-coordinate and Y as the velocity field log-coordinate. 



3.3. The choice of the reference field. To obtain the composition law for the scaling exponents, we 
must make precise the link between A', Y and U, V. This depends on the choice of <f > o and 4>o- As emphasized 
in [1], a convenient choice is to take <f>o =< ln(</>f) > and V’o =< In(V’r) >• Physically, these quantities are 
the most probable values of the fields <j> and 4k The log coordinate X and Y then characterize deviations 
with respect to these most probable values, and are a measure of the intermittency of the fields. In such 
simple case, it is easy to check [6] that the scaling exponents for cj> and 4’ are simply connected to X and Y, 
via: 


0»£(n) -$,£(«) |»=o = ^r, 

dY 

( 3 . 8 ) d P a P )-d P ap)\ P = 0 = ^f. 

Now, the connection between U,V and £(n),C(p) can be made more clearly by exploiting the remark 
made in Section III. A. Indeed, considering the transformation from A = 0,P' = 0 to X n = (<9 n £(n) — 
d n £(n)\n=o)T,Y p = ( d P Cip ) - d p ((p)\ p=0 )T , we get: 


U(n) = cU(n) - cU(n)U =0 = 

dY 

(3.9) V(p) = d p ( (p) - d P C(p)| P =o = 

This shows the existence of a bijection between (U,V) and ( n,p ). It unambiguously defines an internal 
composition law <§ on (n,p), through a transport of the group structure ®: 


(?7,y)((n,p)®(n>') = (U,V)(n lP )®(U,V)(n',p'), 
(n,p)®(nV) = (ff- 1 ,p- 1 )((17,P)(n,p) 
(3.10) ®(U,V)(n',p')). 
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From this property, we can deduce that the shape of 0 is the same as (3.7), with different parameters, 
characterizing the fixed points in the (n,p) space 1 . We have ( n",p ") = (n,p)®(n',p') with: 

n + n' — ann' — Spp' — n'p + np') 

1 — rinn' — XPP 1 

p + p' — Tirin' — ppp' — r6(n'p + np') 

1 — rjnn' — \pp' 

where the parameters a, p,j, 5, t],t, x characterize the fixed points of 0 and satisfy: 

(3.12) V = TX, -\' + <a<i + 7P = 7 2 +<^ 2 t 

3.4. Computation of the scaling exponents. The correspondence established through eq. (3.10) 
enables to classify the statistics of scale invariant random process, i.e., of the possible shapes for £(n), C(p), 
on universal grounds. Consider eq. (3.10) with («, p) = ( n' ,p ') = (1, 1). This provides a recursive equation 
which allows the computation of ( U , V) for any integer m, and then, by continuation, on any real number, 
namely: 

E(m) = (1, 1)®(1, 1)0...0(1, 1) = (1, 1)H 
(U,V)(E(m)) = (C7(1),V(1))0...0(£7(1),V(1)) 

(3-13) = (£7(1),V(1))K 

Here the notation [m] (resp. [m]), stands for the m th iterate via 0 (resp. 0): we have introduced the 
notation E(m) = (E\{m) , E^im)) to generalize the exponentiation m to real, defined by using infinitesimal 
n and p. It is now technically possible to calculate the values of the £(n) and £(p). The set of equations 
(3.7, 3.10, 3.11, 3.13) is a mathematically well-posed problem. Indeed, inverting (3.8), we get: 

an an 

(3.14) ^) = <M|„ 0 + ,, (1) K'W). 

an an 

This provides a symbolic representation of the possible shape for the scaling exponents, depending on the 
fixed points of the composition law. As an example, we treat below some tractable examples, where the 
fixed points are simple. Since we are interested in the case where no coupling exist from the field </> onto 4’ 
(passive scalar limit), we see that this imposes c = t = 0. These values are adopted from now on. 

4. Passive scalar advected by a Gaussian field. We consider here the case where u is a Gaussian 
velocity field, with 2n-th order structure functions scaling like: 

(4.1) (^{Su()‘ 2n '^ = (^{u(x + £) — u(x)) 2 ^ oc 
Because u is Gaussian, it follows the simple “normal” scaling property: 

(4.2) C(2n) = f»C(2). 

We are interested in the scaling properties of a passive scalar, advected by such a velocity field. Specifically, 
we are looking for the scaling exponents £(2n) of the 2 n th order structure function of temperature increments 
defined as: 

(4.3) < (6(x + £)- 0(x)) 2n >~ £^ 2n \ 

Physically, these fixed points correspond to the limit beyond which the moments of the distribution in <f> or i/i are divergent, 
see [6]. 
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In the Kraichnan model [10] of the eq.(l.l) where u is delta correlated in time, a number of exact results 
have been recently established. Kraichnan [10] and Gawedzki and Kupianen [11] have shown the second 
order scaling exponent £(2) is related to the exponent £(2) of the second order structure function of the 
velocity field by: 

(4.4) £(2) = 2 - C(2). 

For scaling exponents of higher order structure functions, both perturbative and non-perturbative results 
are available. In the limit where £(2) — >• 0, it was shown by Bernard et al. [12] that the scaling exponents 
follow the “anomalous” law: 

(4-5) C(2n) = <(2)- 2n ^ n + ~ 2 1) C(2), C(2) -»• 0, 

where D is the space dimension. In the limit where £(2) — >• 2, it was shown by [13, 14] that the anomalous 
correction vanish, so that: 

(4.6) £(2 n) = n£(2), <(2) -»• 2. 

Chertkov et al. [15] used an expansion in the inverse 1/D of the space dimensionality to obtain: 

(4.7) £(2n) = n£(2) - 2n( ^~ 1) C ( 2 ), D _1 -> 0, 

which is consistent with the result (4.5) obtained by Bernard et al. [12]. 

On the other hand, Kraichnan [10] used the linear ansatz for the molecular-diffusion term and found a 
non-perturbative result: 

(4.8) £(2n) = \y/4nDZ(2) + (D - £(2))* - \{D - £(2)). 

Note that this result is inconsistent with the £(2) expansion of Bernard et al [12]: using (4.4), one finds that 
in the limit £(2) ->■ 0, the exponents £(2 n) - n£( 2) given by (4.8) tend to a finite limit which contradicts 

(4.5) . 

We finally mention the exact non-perturbative calculation of the anomalous correction to fourth scaling 
exponent p 4 = 2£(2) — £(4) obtained by Benzi et al. [16] in the case of a random shell-model for passive 
scalar. This discrete toy model mimics most of the properties of the equation (1.1) [17]. In such model, the 
properties of the scaling exponent are independent of ultraviolet and infrared boundary conditions, except 
in the case where £(2) — > 0 if the diffusive scale is kept fixed. In this last case, diffusive effects destroy the 
inertial range scaling, and p 4 —>■ 0, like in the models of [12] or [15], in the limit £(2) — ¥ 2. Chertkov et al 
also shows p 4 = 4(2 — £(2 ))/D [18]. 

4.1. The scaling exponents. Let us now examine those predictions within the framework of the scale 
invariant theory obtained in Section 3. In the stationary situation and in the non-diffusive limit, the equation 
(1.1) is invariant under the transformation Sh,h'{ A). This shows that, for the range of scale in which the 
non-diffusive limit is meaningful (the so-called inertial range), the structure functions obey power laws (4.3), 
where the scaling exponents £(n) may be a priori computed by using scale invariant theory. 

In the Gaussian case, all the moments are convergent. The convergence of all moments implies that 
the nonzero fixed points of the composition law for p must be — oo or oo. That is possible only if r = r) = 
X = p = 0. Since the scaling exponent of the Gaussian field is linear, that is £(2n) = n(( 2), V (n) has to 
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be identically zero according to (3.8). With c = h = 0 as imposed by the passive scalar requirement and 
b = k = 0 required by the linear scaling exponent, the relevant composition law becomes: 


(4.9) 


and 


U" = U + U' - aUU' 

V" = V + V' 


(4.10) 


n" = n + n' — arm' 
p" = p + p' 


Note that it is compatible with the requirement V (p) = 0. The solution of eqs. (4.9) and (4.10) can be 
found by iteration, as explained in Section III.D. We now introduce the parameter (3 via 


(4.11) 


f3 = 1 - a*7(l). 


The solution U(p) of the equation (4.9) takes two different forms corresponding to the two following cases: 
(1) If [3 ^ 1 and a/0, then 


(4.12) U(n) = - (l-(l-an) A_1 ) 

where A - 1 = In (3/ ln(l - an(l)). 

Using eq.(3.8) and integrating over n, we obtain: 

(4.13) £(2n) = 2nh 0 - (l - (1 - 2cm) A ) , 

where ho is a constant depending on <9 n £(0). This expression is a slight generalization of Kraichnan’s formula. 
Indeed, setting ho = 0, A = 1/2, a = -2 (D - £(2))/£>£(2) and a = —2D£(2)/(D - £(2)) 2 , we find exactly 
(4.8). The ’’log-normal” formula of Chertkov et al (4.7) and Bernard et al (4.5), with £(2 n) being the sum 
of a linear and a parabolic contribution is then obtained from (4.13) by taking the limit either A — > 2, or the 
limit as follows. 

(2) (3 = 1 and a = 0. Then, 


(4.14) £(2n)-n£(2) = 2^n(n-l). 

n(l) 

Note that this limit is not achieved when the parameters are “tuned” to Kraichnan’s value in the limit 
((2) — >- 0, £(2) — > 2. This illustrates in another way the discrepancy arises between Kraichnan’s result and 
Bernard et al.’s asymptotic result. 

Finally, we get two additional interesting possibilities when 
(3) j3 = 1 and a/0 


(4.15) 


U(n) = 


U( 1 ) 


ln(l — cm(l)) 
£(2 n) = 2nho + 

U( 1 ) 


ln(l — an), 


aln(l — an(l)) 


[2an — (1 — 2an) ln(l — 2an)j. 
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This expression tends towards the log-normal asymptotic form when a — >■ 0. 

(4) ^ 1 and a = 0 

U(n) = ^(l-/3 n/nW ), 

(4.16) £(2n) = 2nh 0 + (l - /l 2 "/"! 1 )) . 

This is the log-Poisson law [19, 20, 21], which again tends to a log-normal in the limit (3 — ¥ 1. 

4.2. Discussion. Our result obviously generalizes the results of [10, 12, 15], using four arbitrary con- 
stants. Because we use symmetry arguments, we cannot obtain any quantitative information about the 
value of the anomalous correction. Only “matching” with exact computations (starting from the original 
equations) can provide values for these constants. In the present case, we have four arbitrary parameters: a, 
17(1), a and ho- These parameters can be arbitrary function of ((2) which need to be determined. The exact 
result (4.4) enables to write an equation linking our four parameters to C(2). Besides that, the asymptotic 
expressions of Bernard et al. and Chertkov et a 1. only provide asymptotic expressions for 17(1), and the 
condition that a must tend to zero as D — > oo or £(2) — 0. This obviously does not constrain our parameters 
enough, so we shall not bother to try to exhibit one out of the many solutions which fits all the requirement. 
It would be interesting, as one gets more and more precise numerical result as those obtained by Frisch et 
al. [22], to try to derive the four parameters for any given value of ((2), and see if they are enough to fit the 
whole collection of exponents at this value of C(2). This would provide a test of the present theory. 

5. Passive scalar advected by turbulent flow. In this Section, we consider the case where the 
velocity field u follows the Navier-Stokes equations: 

VP 

(5.1) d t u + (u • V)u = bt'Av + f. 

P 

Here p is fluid density, f force, P pressure and v viscosity. Eq. (1.1) indicates that the passive scalar is 
transported by the turbulence. Eq. (5.1) shows that the dynamics of the velocity field is not influenced by 
the passive scalars. 

In the stationary unforced inviscid limit (is, n -¥ 0), the set of equations (1.1, 5.1) is invariant under 
the transformation Sh.h'W- This shows that for the range of scale in which the inviscid unforced limit is 
meaningful (the so-called inertial range) the structure functions of velocity and passive scalar obey power 
laws (1.2). Dimensional considerations lead to the prediction £(n) = C( n ) = n/ 3. However, the experiments 
and numerical simulations do not support the predictions. These exponents are observed to deviate from 
the linear laws and the deviation is higher for the temperature than for the velocity. It is believed that the 
deviations are due to intermittency of velocity and scalar dissipations. 

Many models have been proposed to reproduce the scaling exponents of velocity structures. Among 
them, She and Leveque’s model [19] gives an excellent prediction of experimental results without any free 
parameters. Later, this model is shown to be equivalent to assume that velocity increments follow a log- 
Poisson statistics [20, 23]. The log-Poisson statistics converges to a log-normal statistics under a suitable 
limit. In fact, a recent careful analysis [24] of turbulent data using wavelet transform indicated that the 
turbulence is so close to log-normal that distinguishing between log-Poisson and log-normal statistics is 
impossible, given the available amount of turbulent data. Therefore, we shall assume in the following 
discussion that the turbulent velocity is log-Poisson, and then take the suitable limit towards log-normal 
statistics to examine whether the passive scalar statistics can help in making the distinction. 
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As an extension of log-Poisson statistics of velocity fields, Cao and Chen [25] assumed a bivariable log- 
Poisson distribution for joint moments of velocity dissipation e( and scalar dissipation N{. This assumption 
implies that passive scalars also obey log-Poisson statistics, which is in remarkable agreement with the results 
of numerical simulations. Ruiz, Baudet and Ciliberto’s recent experiments [26] also provide strong evidence 
to support log-Poisson statistics of passive scalars. He et a 1. [21] propose a hierarchical relation for the joint 
moments to explain why passive scalars must be log-Poisson. In the following, we explore the possibility of 
log-Poisson statistics of passive scalars within the formalism developed in Section 3. 

5.1. The scaling exponents. A remarkable simplification occurs from dynamical considerations: since 
velocity dynamics is not influenced by the scalar dynamics, that is, the composition law for V must be 
independent of U, one must choose c=h = T = i) = 0. Moreover, to select the log-Poisson statistics for 
velocity increments [4] requires k = 0 and f3 = k = 0. Therefore, the relevant composition law in this case is: 


(5.2) 


and 


U" = U' + U- aU'U - dV'V - j[U'V + UV'] 
V" = V' + V- bV'V 


(5.3) 


n" = n' + n - an'n - Sp'p — y[n'p + np'} 
p" = p' + p 


A rescaling of U and V by an arbitrary factor t , and of a, b, d and j by l/t leaves (5.2) unchanged. At 
the same time, a rescaling of n and p by l/t , and a rescaling of a, <5, 7 by a factor t leaves (5.3) unchanged. 
Since the scaling exponents depends only on the derivative of U and V with respect to n and p, we may 
assume without loss of generality that, say, p( 1 ) = 1 . 

The equation (5.2) has one stable fixed point 


< 5 - 4 ' u * - ^ * - \ 
and one unstable fixed point 

< 5 - 5) u - = v • = i 

We solve (5.2) and (5.3) and get an implicit representation of the solution via, 


U[n(m )] = (U + - I/_)(l - (1?) + W(1 - P?) 
V\p(m )] = K(l-/?2 ro ) 


where 


(5.7) 

and 


A = i 
A = i 


^(D 

u+ - U- 


V'(l) 

v; ' 
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(5.8) | 

f n(m ) 
l P(rn) 

( a 2 \ 1 -af a 2 m 

= ( n(1) l-Jl-,,*!-,, 

= m 

where 


ci\ = 1 — cm(l) — 7, 

(5.9) 


0,2 = —5 — 7 n(l). 


Here, n( 1) is a free parameters. To avoid any discontinuity while ensuring n(0) = 0, we shall only assume 
that n(0) tends to 02 when 01 tends to 0. Note here that the log-normal shape for the velocity increments 
is obtained by letting 0 2 — > 1, in which case the function U simply becomes: 


(5.10) 


U(m) = (U+ - U-) (1 - 0?) + m. 


To get U(n) as a closed form, one needs to invert the function n(m) given by (5.8). This is analytically 
possible only for special values of the parameters, otherwise one can only use the implicit representation 
(5.6) and (5.8). The special values where an explicit formulation of U(n) is possible are: 

1) when o 2 = 0: we get: 

U(n) = (U+ - U-) (1 - (1 + n/as)^- 1 ) 

+ K(1-(1 + n/a 3 ) A2 " 1 ), 

V{p) = V. (1 - 0 P 2 ) , 

a 3 (U+ - U.) 


(5.11) 


£(n) = nh $ + 
K 03 

C ip) ~ph ~ 


Ai 


(1 - (1 + n/a 3 ) Xl ) 


(1 - (1 + n/a 3 ) X2 : ), 
V, 


In 02 


(1-0 2 P ), 


where A* — 1 = ln/3,/ln(oi), for i = 1, 2, a 3 = n(l)/(oi - 1) and h $ and h ^ are integration constants. The 
shape for £ is the usual log-Poisson shape. The shape for £ is a superposition of log-levy function. In the 
log-normal limit for the velocity increments, one gets a different shape: 

03 (U+ - U.) 


£(n) = nhs + 


Ai 


(l — (1 + n/a 3 ) Xl ) 


(5.12) 


+ o 3 ((1 + n/a 3 ) ln(l + n/a 3 ) - n/a 3 ) , 

C (P) =ph c + y. 


So the log-normal limit is not really simpler than the general solution. 
2) when on = 0, we get: 


U{n) = (17+ - U- 

V(p) = K (i - ft ) , 

n(l)(U + - U-) 


) (1 - /3r /n(l) ) + v * (1 - 0 2 ) n/n{i) ) , 


(5.13) 


£(n) = n/i^ 

_ n(l)K 
In 02 

C(p) =ph C - 


In 0i 

(] V~0 n 2 /n[X) ), 
(l-0 P 2 ), 


(!-/?”/”(!)) 


V, 


In 02 
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where h q and are integration constants. In this case, they are scaling exponents connected with the 
properties of the rarest events (the most intermittent structures), whose connexion with finite size effects is 
discussed in [5]. The physical interpretation of the parameters appearing in (5.13) is obtained by taking the 
limit n — >• oo or p — >• oo, when pi and fc are less than one. We find that 

£(n) ~ n/i£ - n(l)(U+ - 17_)/Infr - n(l)V./ln& 


and 


C(p) <xph c - 14/ In ^2- 


This suggests interpreting — n(l)(U + — U-)/\npi — n( 1)14/ In ^ and — V*/\np 2 as the codimensions of 
the most intermittent structures of the passive scalar and of the velocity. At this stage, we can make an 
interesting remark: when n(l) is small (this is a case of weak coupling between passive scalar and velocity), 
the last term in (5.13c) can be neglected. In this case, the scaling exponents of the passive scalar have a log- 
Poisson shape. We show later by fitting the experimental data that this coupling coefficient is indeed small. 
This therefore explains the observation in [26]. Also, we note that the scaling exponents (5.13) reassemble a 
log-Poisson distribution with two “atoms” , whose existence was conjectured in [23] . Finally, it is interesting 
to consider the limit where the velocity increments are log-normal. In that case, we get: 


(5.14) 


£(n) = n/i$ - 

C (p) = Ph + 


n(l)(U+-U.) 

In Pi 
V{l)p 2 


(i - ft 


n/n{l)\ 


+ 


V(l)n 2 
2n(l) ’ 


In that case, the scaling exponents are the superposition of an exponential log-Poisson shape, and of a 
parabolic log-normal shape. At large positive n (if pi < 1), the parabolic shape dominates the scaling 
behavior, and one does not observe the simple linear asymptotics existing in the log-Poisson case. This 
observation could be used as a tool to discriminate between log-Poisson and log-normal velocity increments, 
by considering the statistics of passive scalars. This is an indirect check, but in some experimental settings, 
this might actually prove easier than the direct check. 

In all other cases, we may only obtain an implicit representation of the scaling exponent by using the 
chain rule d£/dm = d£/dn ■ dn/dm. 


5.2. Value of the constants. We obtained different explicit and implicit forms of scaling exponents 
depending on several constants. In the most general (implicit) case, we have 9 free parameters. In the 
simplest case (number 2), we have 6 free parameters. We need the same number of constraints to compute 
these parameters. For h^, Pi and V* / In p 2 , we can use the values obtained by [19], which provide a good 
fit of the velocity increments: h ^ = 1/9, p 2 = 0.88 and V*/ln/3 2 = 2. For the remaining parameters, 
we need to use values measured for the passive scalar. Given the available precision of the measures, it 
is not really realistic to consider more than n = 5 or n = 6 scaling exponents for computing the values 
of the constants. It is then only possible to determine the constants in the log-Poisson case, where only 
three remaining parameters need to be determined. This simplification enables a partial determination 
of the parameters. If the most intermittent structure have respectively vortex filaments and temperature 
sheets, as observed in numerical simulations [27, 28], their codimensions are approximatively one and two. 
By considering that the most intense energy dissipation events are only driven by the mean transfer term 
(< (5v) 3 >) and using the K62 refined similarity hypothesis [29], She and Leveque [19] obtain h ^ = 1/9, and 
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n 


Fig. 5.1. Comparison between the experimental data of Ruiz et al. (diamonds) and formula (49) with: A set of parameter 
derived by physical arguments copying the argumentation of She and Levegue [19] for the velocity increments. This gives 
(/i£, /3i,n(l), U+ — U- ,h^,Vi,, fo) = (1/9,0.985,0.1,0.015, 1/9,0.26,0.88) (dash line); A set of parameters derived by a best fit 
to the data (hf, /h, n(l), U+ - U-,h 0 K,/3 2 ) = (0.06,0.915,0.25.0.12,1/9,0.26,0.88) ( solid line). 

then, using £(3) = 1, /3| = 2/3. A same kind of argument also leads to h £ = 1/9. Indeed, if we suppose that 
the most intense temperature dissipation events are only driven by < (8v) 3 > and < 5v(59) 2 >, we obtain 
N°° ~ 0 2 < (5v) 3 > 1 / 3 / < Sv(8Q ) 2 >, where 0 is some constant temperature. Since < (8v) 3 >~ l and 
< Sv(80) 2 >~ £, we obtain N 00 ~ l~ 2 ! 3 . Then, by using the refined Obukhov-Corrsin hypothesis ([30]), we 
get \ — 1/9 + hjv/2 + 1/3 = 1/9. We may then use these values and vary n(l) and (i\ to obtain the best fit 
with the data. The result is given in Figure 1 for n(l) = 0.1 and (i\ = 0.86. It is not really good, especially 
at large n where the formula gives a much larger value than the data. This is because the value = 1/9 is 
too high. A three parameter fit of the data gives h $ = 0.06, n(l) = 0.25, (3\ = 0.7, resulting in a codimension 
of 0.84 for the passive scalar rarest events. The result is displayed in Figure 1. 

Certainly, more work is needed to understand the meaning of the parameters. A recent work of Dubrulle 
[5] [31] suggest that some of them might be linked with finite size effects. Work is in progress in that 
direction. 

6. Conclusion and discussion. In this paper, we develop a formalism leading to the equations fol- 
lowed by the scaling exponents of two-dimensional scale invariant random fields. It is the extension of [1] on 
one-dimensional random fields. We do not make any assumption on (statistical) relations of the two random 
fields. For example, they may be independent of each other or anisotropic. Thus, the two random fields pos- 
sess independent scaling exponents and they may have different limit velocities. This situation is completely 
different from the relativity theory of motion, where the limit velocity is the same one (light velocity) in 
each direction. We begin with the structures of transformation group, and then compute the transformation 
matrix. The composition law, derived from the similarity transformation, indicates the equations relating 
the different scaling exponents. Based on the scaling exponent equations, the classification of all exponents 
is possible by the computation of the fixed points of the transformation. This is performed in a simple case, 
the passive scalar advected by i) a random Gaussian velocity field; and ii) a turbulent velocity field. 
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In the Gaussian case, we show that the passive scalar increments follow a log-levy distribution gener- 
alizing the solution of Kraichnan. This distribution tends in an appropriate asymptotic limit, towards a 
log-normal distribution. The anomalous correction to the scaling exponent can be matched with the Kraich- 
nan solution [10] or with recent perturbative expansions of Chertkov et a 1. [15] and Bernard et a 1. [12] via 
an appropriate tuning of the parameters, which cannot be predicted by our simple symmetry arguments. 

In the turbulent case, we show that when the velocity increments follow log-Poisson statistics, the passive 
scalar increments follow various different statistics, the simplest one being close to log-Poisson. It is in fact a 
log-Poisson distribution with two “atoms” , thereby providing an explicit example of a random system with 
more than one atom whose existence was conjectured in [23]. This result provides an explanation of the 
experimental findings of [26] about the temperature increments. 

We note that scale invariance may play a role in many other physical problems. It would then be 
interesting to consider other values of the parameters to describe other physical situations such as scaling 
laws in magneto-hydrodynamics. 

Acknowledgments G. W. He expresses his gratitude to Lachieze-Rey for discussion around the relativity 
of motion. He also would like to thank SAP for its hospitality. B. Dubrulle acknowledges useful discussions 
with M. Vergassola. 

Appendix A: Derivation of Similarity Transformations 


We show that the similarity transformation compatible with the following postulates are the solution 
(3.3) by generalizing the method [1] to the two-dimensional case: 

Posulate 1. The log-coordinates are invariant under global translation. 

Posulate 2: Among all imaginable reference fields, there exist a continuous class of equivalent reference 
fields built on scale invariant fields. 

Posulate 3: the group of similarity transformations is commutative. 

Postulate 2 implies that the composition of two similarity transformations with respective parameters 
(U, V) and (U',V') should be a similarity transformation, associated to a third parameter noted ( U",V "). 
Composition of two transformations in (3.2) then gives 


A" 

ii 

+ BD' + , G'UU' 

B" 

ii 

to 

' + BE' + , H'UV' 

, "U" 

= A, ' 

U' + B,'V' + „ 'U 

D" 

II 

b 

’ + ED' + , G'VU' 

E" 

= DB 

' + EE' + , H'VV' 

, "V" 

= D, ' 

U' + E, 'V' + „ 'V 

G"U" 

= GA! 

'U + HD'V + , G'U' 

H"V" 

= GB 

'U + HE'V + , H'V' 

n 

= G, ' 

UU' + H, 'VV + „ 


We now use the Postulate 3 about commutativity. We first apply the transformation ({/', V') and then 
the transformation (U,V). By comparison with (6.1), we obtain several equations. For example, equations 
(6.1.1) and (6.1.5) give 

BD'-B'D = UU'(,'G-,G'), 

DB'-D'B = VV'{,'H-,H'). 
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Factorizing (6.2) leads to 
(6.3) 


D = cB , G = —h, , H = —k, , 
where c, h , fc are adjustable parameters. 

Substituting (6.3) into (6.1.2) and (6.1.4), and then comparing their coefficients, we obtain 


h = ck. 


The same technique applied to (6.1.3) and (6.1.6) gives 

A = (1-aU- jV ), 

(6.4) B = -( jU + dV ), 

E = (1 - cdU - bV), 

where a,b,d,j are adjustable parameters. (6.3) and (6.4) are solutions of equations (6.1) provided the 
parameters follow 

—k + ad + jb = j 2 + (Pc. 

This completes our proof. 
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